Gigahertz quantized charge pumping in graphene quantum dots 
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Single electron pumps are set to revolutionize 
electrical metrology by enabling the ampere to be 
re-defined in terms of the elementary charge of an 
electron^. Pumps based on lithographically-fixed 
tunnel barriers in mesoscopic metallic systems^j] 
and normal/superconducting hybrid turnstiles^f 
can reach very small error rates, but only at MHz 
pumping speeds corresponding to small currents 
of the order 1 pA. Tunable barrier pumps in semi- 
conductor structures have been operated at GHz 
frequencies^^, but the theoretical treatment of 
the error rate is more complex and only approxi- 
mate predictions are available 8 . Here, we present 
a monolithic, fixed barrier single electron pump 
made entirely from graphene. We demonstrate 
pump operation at frequencies up to 1.4 GHz, 
and predict the error rate to be as low as 0.01 
parts per million at 90 MHz. Combined with the 
record-high accuracy of the quantum Hall effect 9 
and proximity induced Josephson junctionals, ac- 
curate quantized current generation brings an all- 
graphene closure of the quantum metrological tri- 
angle within reaclfGQ Envisaged applications for 
graphene charge pumps outside quantum metrol- 
ogy include single photon generation via electron- 
hole recombination in electrostatically doped bi- 
layer graphene reservoirs 12 , and for read-out of 
spin-based graphene qubits in quantum informa- 
tion processing^. 

The graphene double quantum dot pump used in this 
work consists of two lithographically defined graphene 
islands which are coupled to each other and to source 
and drain contacts by narrow constrictions [(Fig. Qa)]. 
The strong Coulomb interaction between electrons blocks 
the continuous flow of current through the device until 
nearby plunger gates lower the energy cost of tunneling 
through both quantum dots. To pump electrons from 
source to drain we rapidly modulate the voltage on the 
gates such that only a single electron can be transfered 
through the structure per modulation cycled A single 
cycle comprises the three stages which are illustrated in 
Fig. [I] (i) the electrostatic potential on the first quan- 
tum dot is lowered as the voltage Vgi on plunger gate 1 
increases, pulling an electron in from the source lead; (ii) 
the potential on the second quantum dot is lowered as the 
voltage Vg2 on plunger gate 2 increases, shifting the elec- 
tron to the second dot; (iii) the electron is pushed out to 
the drain contact and the system returns to its original 
charge configuration ready to pump the next electron. 




FIG. 1: Adiabatic pumping mechanism in a graphene 
double quantum dot. (a) Atomic force microgaph of the 
device and experimental setup used to measure and generate 
the pumped current. An oscillating voltage Vkf(£) is added to 
the DC voltages Vgi and Vg2, which control the total number 
of electrons on each quantum dot. The phase delay added 
between the modulation voltage applied to Vgi and Vg2 con- 
trols the shape of the pump loop in the space of the two gate 
voltages, (b) DC current through the device as a function of 
DC voltages applied to the gates. Red dashed lines indicate 
the edges of the honeycomb stability diagram (See Supple- 
mentary Section £ |I h\ . A quantized current is pumped when 
the gate voltage modulation produces a trajectory (yellow) 
that encircles a triple point, passing through the sequence of 
transitions (i)— »(h)— >-(iii). 



The frequency / of the oscillating voltage Vrf applied 
to the gates determines the rate at which electrons are 
transferred, and thus the size of the pumped current, 
I = ef. 

Figure. [2]ja) shows the pumped current generated by 
applying a modulating voltage at / ~ 1.5 GHz while 
sweeping the voltages Vgi and Vgi in the vicinity of a 
pair of triple points (Supplementary Fig. S^]). We ob- 
serve extended plateaus in the pumped current at equal 
and opposite values around each triple point [Fig. [2jb)], 
demonstrating that the charge pump delivers single elec- 
trons from source to drain per cycle. The opposite sign 
of the pumped current when encircling different triple 
points in a pair is expected, as the equivalent sequence 
of configurations pumps an electron in the opposite di- 
rection, or holes in the same direction^. Unambiguous 
confirmation of quantized charge pumping is shown in 
Fig. [ljc), which plots the pumped current as a func- 
tion of / with the DC voltages fixed on the right hand 
triple point. Aside from the oscillatory behaviour, which 
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is introduced as the pump loop periodically changes di- 
rection around the triple point (Supplementary Section 
SID I, we observe remarkable adherence to the quantized 
value / = ±e/ over a range of frequencies up to a few 
GHz. At fixed pump frequency and gate voltage we also 
observe a plateau in the pumped current as a function 
of the RF power controlling the size of the pump loop 
[Fig. [2jd)]. The independence of pumping efficiency with 
power confirms that our pump loop can satisfy the purely 
topological requirement of encircling a triple point. It is 
likely that the departure from ef for negative currents 
derives from leakage currents due to intersecting triple 
points in the vicinity of the left triple point [See Supple- 



mentary Section £IA 



Owing to the demand for high currents in the major- 
ity of charge pump applications, a key figure of merit 
is the quantization accuracy at high frequency®. The 
speed of tunneling through potential barriers illustrated 
in Fig. [I] is limited by the RC time constant of the bar- 
riers. To ensure that quantum fluctuations do not de- 
stroy charge quantization on the islands, R 3> Rk- Con- 
ventional metal-oxide junctions formed by electron beam 
lithography have capacitances of order 1 fF and there- 
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FIG. 2: Gigahertz quantized charge pumping. (T=300 
mK) (a) Pumped current generated by the double quantum 
dot as a function of Vgi and Vg2 at / =1.465 GHz (e/ «235 
pA) and P=-15 dBm. (b) Plot showing the pumped cur- 
rent along a line passing through both the triple points, (c) 
Pumped current as a function of frequency (P—-15 dBm) and 
(d) power (/ =1.465 GHz). Red and blue dashed lines show 
the quantized values I — ±e/ of the pumped current. 



fore RC ~ 0.1 ns. In pumps formed from chains of these 
metallic junctions, the error rate due to missed tunnel 
events can be as low as 0.01 parts per million for / < 
10 MHz, but is expected to degrade exponentially with 
increasing frequency, approaching 1 for / = 1 GH d 2 l 17 l 
On the other hand, tunnel barriers formed in silicon®^ 
and graphene^ have capacitances at least an order of 
magnitude less than metal-oxide barriers with the same 
tunnel resistance, leading to faster RC time constants, 
although high speed quantized pumping in a semiconduc- 
tor fixed-barrier pump has not previously been demon- 
strated. In our graphene device we estimate R ~ 100 kil 
from the differential resistance at high bias voltage, and 
C ~ 40 aF from the C oulomb charging energy (Supple- 
mentary Section £ I B I . 



In Fig. [3] we show the pumped current as a function of 
frequency for a different pair of triple points to the data in 
Fig. [2j At low frequency (lower right inset), the current 
is equal to ef, within the ~ 0.2% calibration accuracy 
of our ammeter, but at high frequency the current falls 
below ef, a behaviour we attribute to an increase in the 
number of failed pumping cycles. In the upper left inset 
we plot the error AI/I = (ef — I) /I as a function of cycle 
time t p = 1//, for the higher frequencies where the error 
is substantially larger than the measurement uncertainty. 
A fit of this data to the expected exponential dependence 
AI/I oc exp[—t p /T] yields r ~ 0.6 ns. Similar fits to er- 
ror measurements on 5-j unctiorP^ and 7- junctiorP metal- 
lic pumps yielded t ~ 6 ns. This scaling of the observed 
error rates is expected from the roughly 10 times scaling 
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FIG. 3: Frequency dependence of the pumped current. 

Current as a function of drive frequency (1.2 K); Dashed line 
(red) is the linear relation expected from pumping one elec- 
tron per cycle. The lower right inset shows the low frequency 
pumped current and the upper left inset shows the percent- 
age error in the pumped current as a function of the pump 
cycle period. See Supplementary Section S;I D| for explanation 
of the frequencies where the current is not quantized. A sys- 
tematic offset of 40 fA due to the uncalibrated electrometer 
has been subtracted from the raw data. 
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FIG. 4: Quantization accuracy at low frequency, (a) Map of the pumped current as a function of the DC voltages Vgi 
and Vg2 while modulating Vrf at / =15 MHz. (b) Line profiles of the pumped current as a function of Vgi at / =15 MHz 
(pink, green) and / =20 MHz (black). The red/blue dashed lines are the values expected from pumping one electron/hole per 
cycle, (c) Pumped current as a function of the pump loop size with the gate voltages fixed next to the right-hand triple point. 



of the junction capacitances, and extrapolating to lower 
frequencies, we predict an error rate due to missed tunnel 
events of less than 0.01 ppm at / = 90 MHz. We note 
that the error rate due to missed tunnel events depends 
only weakly on the number of junctions in the pump, in 
contrast to co-tunneling errors which are strongly sup- 
pressed by adding extra junctions^. 

In Fig. H| we examine in more detail the low-frequency 
pumping and its robustness against variations in control 
parameters. The pump map around one pair of triple 
points at f—15 MHz is shown in Fig. |4|a). The phase 
difference between Vgi and Vqi was optimized so as to 
avoid any DC current contribution (see Supplementary 
Section £1 A for more details of this procedure.) Figure 
|4|b) shows comparisons between the line sections I (Vgi) 
passing through each triple point and and the values 
7=±e/ expected if one electron is transfered per cycle. 
Flat plateaus in I (Vgi) persists with an average value of 
3.189 pA and the extent of the plateau is AVgi«6 mV, 
which is close to the average TP spacing of AVpp~ 5.5 
mV [see Supplementary Fig. £j7] ■ Fig. [4^c) demonstrates 
how robust this quantization is against variations in the 
size of the pump loop. At low Vrf the pump loop neither 
encloses or passes through a triple point, nor intersects 
the cotunneling and polarization lines spanning between 
them, so the time-averaged current is close to zero. As 
the pump loop expands it first passes through a TP and 
the current begins to increase. Once beyond, the pump 
loop encircles the triple point and a current plateau de- 
velops and extends AP«3 dBm, or AVrp-«12 mV. At 
high power the pump loop encloses both triple points 
and the current consequently decreases. These plateaus 
suggest the pump is robust against low-level charge fluc- 
tuations occuring over several hours of normal pump op- 
eration, although ocassional instabilities did occur in cer- 
tain ranges of gate voltage (see Supplementary Fig. S^]) . 
Such instabilities in the charge configuration of the dots 



due to charge fluctuations in the Si02 can be minimized 
by using boron nitride^ substrates which are free from 
dangling bonds and charge trappy. 

We ascribe the high performance of our adiabatic 
pumps to a number of factors deriving from graphene's 
unique two-dimensional physical and electronic struc- 
ture. Firstly, the presence of strong edge and poten- 
tial disorder in lithographically defined graphene nanos- 
tructures leads to the formation of multiple quantum 
dots in the constrictions acting as tunnel barriers be- 
tween the dots (Supplementary Section £IC). Rather 



than this impeding pumping, the resulting capacitance 
and resistance of the random tunnel junctions between 
the quantum dots promotes a high intrinisic tunneling 
time while simultaneously suppressi ng co tunneling events 
due to the large overall dissipation^!!. Secondly, the 
large interdot capacitance and linear electronic disper- 
sion in graphene leads to a large and occupancy depen- 
dent single-particle energy spacing near the Dirac point, 
A(N) = hv F /(dy/N) (N < 5P3 which suppresses pho- 
ton assisted interdot transitions^ and protects adiabatic- 
ity of the pump when operated at high frequency. This 
permits the interdot capacitance to be optimized through 
increased dot size, thereby improving the tolerances for 
parallelization, without compromising on accuracy. The 
large separation between triple points in strongly ca- 
pacitively coupled dots also suppresses leakage currents 
while crossing the polarization line of the pumping cycle 
[(ii)— >(iii) in Fig. |TJb)] . In future, the precise origin of 
error mechanisms at high frequency and the role of co- 
tunneling can be te sted by equipping our devices with 
charge detector d 2 l 27 l 

Our predicted performance of 0.01 ppm at «100 MHz 
implies ten pumps operated in parallel would deliver 
100 pA with metrological accuracy. The monolithic 
processing offered by large-ar ea grap hene grown by de- 
composition of silicon carbide^ 28 * 29 ! or chemical vapour 
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deposition^ is ideally suited for such large-scale paral- 
lelization and integration with other graphene nanoelec- 
tronic devices. In particular, interfacing charge pumps 
with edge states in the quantum Hall regime would lead 
to a high accuracy voltage standard operab le at high 
temperature and lower magnetic field d 9 * 31 *^. Since in- 
duced superconductivity and Josephson effecti^ persist 
in graphene at magnetic fields compatible with the quan- 
tum Hall effectP^, a realization of the quantum metrolog- 
ical triangle in a single graphene device is also now within 
sight. 

The availability of graphene charge pumps will improve 
the prospects for graphene in such areas as single electron 
logic and quantum light generation. Pumping single elec- 
trons into gate defined p- and n-type regions of gapped bi- 
layer graphene^ would generate a stream of tuneable sin- 
gle photons via electron-hole recombination^. Because 
of the close analogy between Dirac Fermions in graphene 
and photons, the single-electron charge pumps presented 
here could play a role in electron optics analogous to 
a single photon source in quantum optics. Combined 
with the ballisticity of electrons in suspended graphene 
or graphene on boron nitride^, this would enable test- 
ing optics elements, such as beam splitters and Veselago 
lenses^, in the single-electron regime. Our pumps can 
be tested against quantum pumping concepts which ex- 
ploit graphene's unique bandstructure to achieve more 
exotic single charge and spin manipulation devices^ 3 -^, 
and also pave the way towards the high frequency manip- 
ulation of charge required for developing graphene spin 
qubitsPl. 

This work was financially supported by the European 
GRAND project (ICT/FET, Contract No. 215752) and 
a EPSRC/NPL Joint Postdoctoral Partnership. 



I. SUPPLEMENTARY INFORMATION 

A. Locating triple points suitable for optimised 
pumping 

As discussed in the main text, the excursion (N\, 
AfeJ-K-Ni+l, N 2 )^(Nx+l, N 2 )^(N 1 , N 2 +l) in gate- 
voltage space is all that is topologically required to pump 
a single electron through the double quantum dot (DQD) 
per cycle, where N\ and N 2 are the number of electrons 
on quantum dot 1 (QD1) and quantum dot 2 (QD2), re- 
spectively. However, for accurate quantized pumping it 
is also necessary to avoid regions of the stability diagram 
where there is a background (non-pumped) contribution 
to the total current, i.e., upon crossing the edges of the 
honeycomb during stages (1) and (3), and between the 
triple points (TPs) during stage (2) [see Fig. (c)]. To 
find a suitable range of gate voltages where it is pos- 
sible to satisfy these conditions by tuning the shape of 
the pump loop, prior to pumping we fix the side gate at 
the charge neutrality voltage^HS] and measure the DC 
current as a function of the DC voltages Vgx and Vgi- 



Fig. ^5] illustrates the variability of this DC contribu- 
tion at T=1.2 K, over a ss2 V range of the gate volt- 
ages, corresponding to «30 single-particle states. Wc 
find large variations in both the size of the high current 
regions close to the vertices of the honeycomb cells and 
the low current ridges spanning their edges. In partic- 
ular, at large negative voltages (Fig. ^5] region A) the 
anti-crossings of the honeycomb structure are less pro- 
nounced, reflecting stronger inter-dot tunnel coupling, 
while at lower voltages (region B) the edges of the hon- 
eycomb are not visible, indicating cotunneling is strongly 
suppressed. These features of the stability di agram h ave 
been observed frequently in graphene DQDiPSEUH] and 
are attributed to the changing transparency of resonant 
states in the constrictions, which form in the narrow 
barriers due to potential and edge disorder and quan- 
tum confinement! 25 * 38 * 39 * 42 * ^. Their presence is betrayed 
by the fact that weak and strong coupling regimes lie 
along bands (pink lines) roughly parallel to the DQD 
resonances (green lines). Such bands are naturally inter- 
preted as broad resonances of smaller localized states^, 
which have similar lever arms due to their symmetrical 
position and proximity to the gates. Their impact on 
operating the DQD as a charge pump is to exclude cer- 
tain regions of configuration space where, for instance, 
single-dot behaviour is observed and the time-averaged 
DC current dominates. The pumped current over the 
same 2 V of gate voltages is shown in Fig. ^6] Although 
regions where the barrier states are more opaque consist 
of absent and irregular TPs, the pumped current is well 
quantized owing to the suppressed cotunneling (absence 
of current along edges of the honeycomb), and is achieved 
without too much optimization of the pump loop. A 
TP pair from this region was used in the high frequency 
pumping measurement shown in Fig. [2^c). 

In the regions where the honeycomb is better devel- 
oped the precise trajectory taken in configuration space, 
which is determined by the amplitude Vrf of the drive 
and the phase shift <fi between the two gates, has to be 
finely tuned. We select region C in Fig. ^6] where the 
honeycomb lattice is uniform to illustrate the pumping 
mechanism in Fig. ^7] With the AC drive on the gates 
(/=12 MHz, P=-25 dBm), we observe overlapping el- 
lipses with crescent-shaped areas of positive and negative 
current that develop around the left and right TPs of each 
pair in the honeycomb [Fig. iQd)]. To understand the 
form and structure of these features, Fig. S^e) shows a 
direct comparison of the locations in gate space around 
a typical pair of TPs. The flat regions (P + and P_) 
within the crescents coincide with where the pump loop 
avoids the DC contributions, whereas the surrounding 
regions exhibit a current \I\<ef owing to the additional 
time-averaged current. The central region (P ) where 
J«0 corresponds to encircling both TPs, leading to re- 
peatedly increasing and decreasing the occupancy of each 
QD without any net transfer of electrons from source to 
drain. Changing the relative phase <fi between the gates 
alters the trajectories around the TPs and consequently 
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the DC contribution to the pumped current. 



D. Oscillations in sign of pumped current with 
frequency 



B. Quantum dot characteristics 

In order to characterize the DQD structure we measure 
the well developed honeycomb structure and bias trian- 
gles at a finite bias of 2 mV, and the source-drain gap 
as a function of Vgi in Fig- The maximum source- 
drain bias gap is AVsd ~6 mV. Assuming that the DQD 
can be modelled as a three junction structure, we expect 
the total capacitance of each dot C% to be related to 
the source-drain gap by AVsd=^g/2Cs- From this we 
obtain Cs ~40 aF, which is typical of graphene QDs in 
the 100 nm size ranged The average gate-voltage sep- 
aration between adjacent charge states is AV G1 =48 mV 
(C G1 =3.3 aF) and V G2 =59 mV (C G2 =2.7 aF), confirm- 
ing that the dominant contribution to the capacitance is 
from the source-drain contacts. To test whether these 
measured capacitances are consistent with dots in the 
lithographically designed positions, we performed elec- 
trostatic simulations using COMSOL multiphysics to cal- 
culate the capacitances and found good agreement with 
our observations. 



To encircle a triple point in configuration space we use 
a single RF generator to modulate the voltages on each 
gate [Fig. iQa)]. Owing to the different cryostat wiring 
and circuitry used to control each gate, at a given fre- 
quency there is a built-in phase difference which leads to 
an arbitrary shape of the pump loop. The path length 
difference to each gate and non-identical frequency re- 
sponse of components, such as the bias tee, introduces 
a frequency dependence to this phase difference. Conse- 
quently, as the frequency is swept in Fig. |ljc) , the phase 
periodically walks-through 2-7T every time the frequency 
increases by 50 MHz, leading to a corresponding oscilla- 
tion in the direction of the pump loop and the sign of 
the pumped current. In order to compensate for this and 
to tune the pump loop shape for optimal pumping, we 
introduce an additional time delay to the RF component 
of Vg2 using a programmable delay line [Fig. SQa)]. Fig. 
^9] shows the pumped current as a function of frequency 
and time delay. The observed decrease in the period 
for current reversal with increasing time delay is entirely 
consistent with our description. 



C. Role of disorder 

Strong edge and potential disorder in lithographically 
defined graphene nanostructures leads to the formation 
of multiple quantum dots in the constrictions acting as 
tunnel barriers between the dotsPEUMIIIIIo] Multiple 
tunnel junctions (MTJs) between localized states were 
also observed in conventional metal-island charge pumps, 
and has in the past served two active roles in single elec- 
tron turnstiles and pumps. Firstly, since turnstile de- 
vices only require a Coulomb gap either side of the QD, 
disordered MTJs can be used despite the fact that their 
precise form and distribution is not knowrPS^ll. The sec- 
ond motivation for incorporating MTJs in pumps stems 
from the fact that cotunneling is suppressed in propor- 
tion to the number of barriers between the source and 
drain. For metrological applications this is particularly 
important as unwanted leakage cu rrent d ue to cotunnel- 
ing destroys quantization accurac y 2 * 17 * 21 !. However, be- 
cause such pumps are operated by clocking an electron 
sequentially through adjacent dots, individual gates must 
be designed to couple strongly to a specific island in the 
MTJ, and the randomly located junctions induced by 
disorder are therefore unsuitable. The high performance 
of the R-pump, which achieves the same error rate as 
MTJ pumps by fabricating on-chip resistors inline with 
the pump to dissipate energy and suppress cotunneling^, 
implies that the disorder commonly observed in graphene 
would improve on the basic pump as it creates a more re- 
sistive environment during pumping. 



E. Dirac point location 

In the absence of a native bandgap, charge localiza- 
tion in graphene relies on lateral confinement to create 
tunnel barriers and a combination of potential and edge 



disorder to produce longitudinal confinement. Fig. E 10 
shows the source-drain current measured at T—1.2 K as 
a function of voltage applied to the side gate. The overall 
form of I{Vsg) has the 'V shape expected for graphene, 
but with a region of suppressed conductance around the 
Dirac point, which is characteristic of graphene nanos- 
tructures at low temperature. The fact that the Dirac 
point is close to Vsg ~ 14 V suggests that there is some 
net charged-impurity hole doping. The transport gap 
AVsg ~10 V reflects the presence of large disorder poten- 
tial fluctuations. Indeed, within the transport gap there 
are many irregularly spaced peaks in conductance due to 
the alignment of energy levels in disorder-induced local- 
ized states. The energy scale of the disorder in our device 
can be estimated from the corresponding range of Fermi 
energy (in monolayer graphene) AEp = hvp\/2a/\e\)^ 
where a is the side-gate capacitanceper unit area, and 
v F = 10 6 ms _1 is the Fermi velocitj^. AE F ss 60 meV 
is typical for exfoliated GNRs on Si02 substrates. Addi- 
tional measurements confirming the location of the Dirac 
point are presented in Fig. £ |11[ which shows the side- 
gate voltage dependence of the Coulomb gap and the 
magnetic field dependence of the transport garJ^'. The 
peak of the characteristic envelope of the charging energy 
and the field dependence of the transport gap symmet- 
rically clos es aro und Vsg ~14 V, indicating this is the 
Dirac poinPSM] We found it was possible to pump over 
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a range of w 2V either side of the Dirac point. 



F. Device fabrication 



Our device is fabricated from a graphcnc flake obtained 
by mechanical exfoliation of natural graphite onto a Si 
substrate capped with a 300 nm thick Si0 2 layer. We use 
an undoped Si substrate in order to minimise stray capac- 
itance, which is crucial for operating the pump at high 
frequency. The DQD structure, comprising two ?s200 
nm diameter QDs separated by short «100 nm constric- 
tions, is etched using e-beam lithography and an oxygen 
plasma. The absence of a back-gate means the electro- 
static potential along the DQD is entirely controlled by 
the voltage applied to the graphene side gates: the volt- 
ages Vci and Vq2 on the two plunger gates are designed 
to change the occupancy of the QDs, while the three side 
gates on the opposite side tune the overall Fermi level. 
To operate the device as a charge pump, the output of 
a single RF generator is added to the DC bias of each 
plunger gate. 
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FIG. 5: Zero-bias stability diagram. (T=1.2 K) Outlines shows the regions where (A) pronounced single dot behaviour is 
observed due to the enhanced transparency of the interdot barrier, (B) the barriers are opaque, and (C) double dot behaviour 
dominates. The triple points shown in Fig [TJb) are captured in the vicinity of region (B). Pink lines run parallel to the 
resonances in the barriers which give rise to the different tunnel coupling strengths, and green lines indicate the slope of the 
single-particle states. 




FIG. 6: Quantized current stability diagram. (/=12 MHz, T=1.2 K) Outline shows the region in Fig. £j7|d). 
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FIG. 7: Experimental setup, (a) Atomic force microgaph of the device and experimental setup used to measure and generate 
the pumped current. An RF signal is added to the DC voltages Vgi and Vg2 used to control the potential on QDi and QD2. 
A phase delay <j> is added to the modulating signal applied to plunger gate 2 and controls the trajectory followed in gate space. 

(b) Plot showing the amplitude of the two plunger gate voltages and the three stages where charge rearrangements occur. 

(c) Double quantum dot stability diagram showing the honeycomb cells where charge configurations are stable, (d) Current 
through the DQD in region C of Figs. ^2j^3]without (left) and with (right) an AC drive added to the plunger gates. A cell of 
the honeycomb has been overlaid for clarity, (e) Plot showing a direct comparison between the DC and AC current behaviour. 
The current polarity and quantization accuracy is determined by the trajectory followed by the pump loop around the TPs. 
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FIG. 8: Stability diagram at finite bias (upper panel) and Coulomb diamonds (lower panel). Upper green line 
shows the gate voltage where the lower panel was measured and the lower green line shows the source-drain bias voltage at 
which the upper panel was measured. Pink lines indicate the separation in gate-voltage between consecutive charge states. 
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FIG. 9: Pumped current as a function of frequency and the added time delay between the oscillatory components 

of Vai and Vg2- 
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FIG. 10: Current as a function of back-gate voltage showing the suppression of conductance around the Dirac 
point with and without an out-of-plane magnetic field (T =1.2 K). 




FIG. 11: Source-drain bias and magnetic field dependence of transport gap (T =1.2 K). 
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